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We investigate the use of an operatorial basis in a self-consistent theory of large amplitude col- 
lective motion. For the example of the pairing-plus-quadrupole model, which has been studied 
previously at equilibrium, we show that a small set of carefully chosen state-dependent basis opera- 
tors is sufficient to approximate the exact solution of the problem accurately. This approximation is 
used to study the interplay of quadrupole and pairing degrees of freedom along the collective path 
for realistic examples. We show how this leads to a viable calculational scheme for studying nuclear 
structure, and discuss the surprising role of pairing collapse. 
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(D . I. INTRODUCTION 

Q 

There is a general quest for understanding of complicated phenomena in terms of a limited set of degrees of freedom, 
chosen through some method appropriate for the problem at hand. Many approaches are available, in areas ranging 
from field theories to atomic physics (see, e.g., the reviews in Ref. 1]). These are typically based on the concept of 
' "relevant degrees of freedom" , or on the introduction of collective motion and collective paths - which are two ways 
^ . to express rather similar principles! 

More specific to the nuclear problem studied in this paper, the old question "what is the correct choice of collective 
coordinate in a many-body system" has had quite a few partial answers, see the review for a discussion of some 
of these. The holy grail of this approach is a method that determines a collective path self-consistently, based only 
on knowledge of the Hamiltonian governing the system. Preferably the method chosen should allow us to measure 
' whether the limited dynamics in a few coordinates makes sense at all, or in the language used above, address the 
' question "how effective are the effective degrees of freedom?" . 

. The constrained Hartree-Fock-Bogoliubov method is commonly used to describe collective paths in nuclear physics 
"'Y' ' (see e.g. 0,1^]). This approach, where the collective subspace is generated by a small number of one-body constraints, 
also goes by the name of generalised cranking. The one-body constraints usually consist of a few carefully chosen 
multipole (particle-hole) operators, as well as a few generalised pairing (particle-particle) ones. For large scale realistic 
^ \ problems such as the description of nuclear fission the number of generalised cranking operators needed in order to 
1 make a realistic calculation becomes very large. There is also no reason to limit the constraints to the standard 
.5^ ' choices; other degrees of freedom, especially those involving spin-orbit interactions might also be important. A more 
satisfactory method should allow the cranking operators to be determined by the nuclear collective dynamics itself. 
$H One such approach, followed in this paper and set out in detail in the review paper Q (a similar approach, plus 

9^. 1 relevant references, can be found in Ref. (Sj), leads to a very well-defined approach, which can in principle be solved 
knowing the Hamiltonian and model space. To find the adiabatic collective path we use the local harmonic approxi- 
mation (LHA). It consists of a constrained mean-field problem that needs to be solved together with a local random 
phase approximation (RPA), which determines the constraining operator. This approach lacks practicality, since the 
size of the RPA problem is, for a system with pairing, proportional to the size of the single-particle space squared. 
Even though enormous matrices can routinely be diagonalised on modern computer systems, the algorithm requires 
repeated diagonalisation of such a matrix, which makes an implementation in realistic calculations prohibitively time 
consuming. 

This requires a solution, or at least a good approximation, and this is the subject of the present paper. A first 
approach to solving the problem has been suggested in the work of one of the authors The best way to test 
such ideas is to use a semi- realistic model, where approximations can be tested against the full method, such as the 
pairing-plus-quadrupole model as employed by Barranger and Kumar in their seminal work. It has been shown 
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in the past |g that at equihbrium the RPA can be solved quite efficiently using a simple basis of operators. Related 
work by Nesterenko et al. |^ may also have some bearing on this problem, but will not be investigated here. In 
short, the idea is that the basic operators of the model, weighted by a suitable power of the quasi-particle energies, 
give excellent results. The state dependence induced by the quasi-particle energies is crucial to the success of the 
approach, and is the main difference with methods based on "naive" constraints. Since the original work was only 
done at equilibrium, we must still check that such a basis of operators provides a good solution along the collective 
path, and we indeed find some important modifications to the method discussed in Once a collective path has 
been found we can diagonalise the collective one-body Hamiltonian along this path, including all the zero modes 
arising from broken symmetries. This will give information on how the collective motion influences the ground state 
properties of the nuclei. In several of the examples discussed below we find low lying states with pairing collapse 
which influence the collective behaviour of the system. If we now quantise the collective dynamics, we must include 
the pairing rotations. This is due to the fact that a point with collapsed pairing behaves similarly to the origin in 
polar coordinates, with the pairing phase playing the role of the polar angle. In this work we shall only study the 
effect of the pairing-rotational modes, ignoring collective rotation for the time being. 

The paper is organised as follows. In Sec. we briefly present the basic principles of our approach, to highlight 
those issues that will make the results easier to understand. The practical form of the equations for the type of 
many-body problem considered here is also discussed, and the form of the approximation is introduced. Results are 
then given in Sec. lIIII and finally we draw some conclusions in Sec. II VI 

II. FORMALISM 

The formalism, as set out in detail in , is based on time-dependent mean field theory, and the fact that a classical 
dynamics can be associated with it. The issue of selecting collective coordinates, and determining their coupling to 
other degrees of freedom, thus becomes an exercise in classical mechanics. Furthermore, if we assume adiabaticity, a 
slow motion where the Hamiltonian can be expanded to second order in momenta, we have a problem that can be 
solved. The solution can be stated without any direct reference to the original nuclear many-body problem and the 
choice of the interaction. 



A. Local harmonic approximation for the collective path in the adiabatic limit 

We assume a classical Hamiltonian depending on a set of real canonical coordinates, ^" (a = 1,...,7V), and 
conjugate momenta, ixp {j3 — \, . . . , N), of the form and vr thus parametrise j^*) 0) 

n{£„^) = {MH\^). (1) 

We shall use a tensor notation, where we use upper indices for coordinates and lower indices for momenta. When the 
same symbol appears as both upper and lower index there is an implicit sum over that index. 

The potential V{C) — H{tt = 0) and the mass matrix B°'^ are given by an expansion of 7i(^,7r) in powers of tt in 
zeroth and second order, respectively, 

'H{^,7r)^V{0 + lB'^^7ra7rp + O{n^). (2) 

Terms of higher order (such as tt^) are supposed to be negligible. The kinetic energy in the Lagrangian formalism 
contains the inverse Bap of the mass matrix, K — ^£,'^Bai3£,^ , and can be interpreted as an inner-product in the 
tangent space to a curved manifold. The inverse of the mass matrix i?"^, is thus the metric tensor; in other words 
the matrix Bap represents the Riemanian geometry in configuration space, since it measures lengths in the tangent 
space. This clearly would not be the case if we had higher order terms in the kinetic energy. 

The central part in our approach to large amplitude motion is a search for collective (and non-collective) coordinates 
which are obtained by an invertible point transformation of the original coordinates preserving the quadratic 
truncation of the momentum dependence of the Hamiltonian |17| , by 

q^ = F{0, {fi,a^l,...,n), (3) 

and the corresponding transformation relations for the momenta and tTq,, 
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where we use a standard notation for the derivatives, 5" = -^g" and f^^ = -^f^- The adiabatic Hamiltonian, 
Eq. ||2Jl, is then transformed into 

7Y(g,p) = V{q) + ^B^^p^p. + (5) 

in the new coordinates. The new coordinates are now to be divided into three categories: the coUective coordinate 
q^, the zero-mode coordinates , I — 2, . . . , M + 1, which describe motions that do not change the energies and 
finally the non-collective coordinates q°' , a = M + 2, . . . ,n. [The approach can easily be generalised to include more 
than one collective coordinate, but that will not be discussed here.] 

The collective coordinate is determined by means of the solution to the local harmonic approach, which consists of 
a set of self-consistent equations. These are: 

1. The force equations 

W,a=A/,„ + Az//„, (6) 

where are the zero-modes (also called Nambu-Goldstone or spurious modes) and A/ represents a set of La- 
grange multipliers (which in nuclear physics are usually called cranking parameters). A is a Lagrange multiplier 
for the collective mode, stabilising the system away from equilibrium (we shall often denote it as the generalised 
cranking parameter) . 

2. The local RPA equation 

VcjB-'f'f^p^ihnf f,o., (7) 

where the covariant derivative V-af3 is defined in the usual way {Vaf3 — (Ka),/9)i 

V^p ^ V^p-TlpV„ (8) 
F^^ = B'<' {Bsp,^+Bs^^p-B^''''). (9) 

Zero modes correspond to zero eigenvalues of the RPA. In principle great care needs to be taken to have 
zero modes behave correctly away from equilibrium. The symplectic RPA j2| is the correct way to do so; 
unfortunately it is rather cumbersome, and as a practical approximation we shall ignore the corrections arising 
from this approach here. In this paper we will also neglect the covariant corrections to the RPA, since they are 
time-consuming to calculate. This means that we do not treat the zero-modes absolutely correctly. 

The collective path is found by solving Eqs. © and lO self-consistently, i.e., we look for a path consisting of a series 
of points where the lowest non-spurious eigenvector of the local RPA equations also fulfils the force condition. In the 
minimum of the potential the spurious solutions decouples from the other collective and non-collective solutions. When 
we are following the collective path away from the minimum one can use special techniques, called the symplectic 
version of the theory, to avoid mixing of the spurious solution and the physical solutions [^] . This has the disadvantage 
that it is numerically much more difficult to implement. In this paper we have chosen to ignore the effects of the 
spurious admixtures to the RPA wave- functions, but these are expected to be small, at least close to the minimum. 
As a result there will be a finite overlap between the collective coordinate and the spurious operators. The price paid 
for these approximations is that at points where RPA frequencies should cross we get narrow avoided crossings. The 
narrowness is a measure of the severity of the truncation errors. One way to by-pass such problems, is to use a basis 
of operators, where such crossings are extremely rare. 

B. Large amplitude collective motion with local harmonic approximation 

The local harmonic approximation has been described in 0]. There the structure is discussed in great detail, as 
is the transition between nuclear physics and classical mechanics. Here the formalism is extended to include pairing 
and constraints on particle number. We start with the time-dependent Hartree-Fock-Bogoliubov equations; in this 
case one finds that a natural choice for the coordinates ^ and tt are the real and imaginary parts of the generalised 
density matrix in its locally diagonal form where the change in the pairing density, tCqq' = ($| aq'ttq |$), can be 
parametrised as 



(10) 
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We want to find a solution of the local RPA equation, J"®™, at the generalised density T?,'^'''^ satisfying the generalised 
cranking equation © 

r — n,p 

where /"°* is also an eigenvector of the RPA equation at 7^"°^ and A'n^p are the particle number operators for 
neutrons and protons, respectively. We use / at the previous point as input, and try to find a point a fixed length 
AQ further along the path, that satisfies Eq. Hll|l for the "old" value of /. Subsequently, a new / is found by solving 
the RPA equations, and this procedure is repeated until Eqs. and l(TT|l are satisfied simultaneously. 
The cranking equation 



,-A("VSr^'- 5: t.N^:i,=Q, (12) 



n 

L J qqi 

r— n,p 

is solved with the additional constraint that 

AQ - (/(") + /(")) • (nl^^ - 7^(")) , (13) 

is fixed (• represent a scalar product). The initial values /^°'' and T?.^*^-' are the results obtained at the previous self- 
consistent point on the collective path. AQ is a measure of the step length in the collective coordinate and Eq. H13|) 
is actually a linear approximation to the integral definition of the change in collective coordinate 

AQ= / Tr(/57^). (14) 

The value of AQ depends on the normalisation of /. In the following we choose the normalisation in such a way that 
the collective mass, B, is position independent, 

B = UB'^^ff^ = 1. (15) 

Equation H12() is solved iteratively by a constrained minimisation, where the change of the generalised density in the 
ith step of the mean-field iteration is given by 

/■(„) Ar20(n) 

- ^i-l'^ + '^i f(n).f(n)^ 20(n) 20(n) + ^-Li y^^) 

■' T=n,p ^^ir ' -"^ir 

T—n,p ) 

where A,^"^7^ = 7^|"^ - n^°\ The step length in the mean- field iteration 5^^^ is chosen to be small for small i, to 
make sure that the iterations converge, but can be chosen larger as the iteration approaches the minimum of the 
constrained mean-field. The parameters e, 77, A and are calculated from the set of conditions discussed below. For 
each n the z-iteration is initiated by choosing 

^0 !<- - ^0 Hn) . Hn) + '^Or 20(n) ^20(n) ' ^^^1 

There are two types of constraints that give the undetermined parameters in the method described above: the fixed 
size of the steps in the collective coordinate and the constraint on particle number. The particle numbers are 
constrained by requiring that AT?, does not change the expectation values N^. Such a constraint can be written in 
diff'erential form as 

^(n)^ ^20(n) ^ ^^g^ 

where t = n,p. We also have to constrain A^'^'t?. in a similar way 

/^i^^-Jl . f(n) ^ 0, (20) 

A|L"^7^ • iVf = 0. (21) 
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The six conditions (Eqs. (|13|l and H19H21|l 'l give a set of equations which can be solved for the six parameters e^""*, 
r{l^\ A^"^ and /i^-"^ for each i and n. The expressions for e, t], A and /i can be found in the appendix. 

The quality of the collective coordinate found above can be quantified and calculated. This is done by calculating 
the decoupling measure, derived in One way of calculating D is by computing the inverse mass matrix Bap 
and then calculate 

where we can approximate the derivative by the finite differences 

d(," /77A7^" , , 

— = V2 . (23) 

The decoupling measure is then calculated to be 

D = Bi^- 1, (24) 

where we have uses the normalisation Eq. H15|) . This quantity is straightforward to evaluate, but it is easier to 
understand from an alternative expression for D, which is based on the generalisation of Eq. H13|) to all coordinates 

Ag" = V2A7^ • (25) 

This leads to 



i.e., D is the sum of squares of the change of the non-collective coordinates with the collective coordinate. This is 
clearly zero for exact decoupling. 



C. Projection basis for the LHA 



One of the main difhculties of applying the LHA method to realistic nuclear problems is the effort required in 
diagonalising the large-dimensional RPA matrix repeatedly within the double iterative process. To limit the compu- 
tational effort we use the method presented in Ref. to reduce the size of the RPA matrix. There it was shown 
that the RPA equation can be solved with good accuracy by assuming that the RPA eigenvectors can be described 
as a linear combination of a small number of state-dependent one-body operators. The quality of the results, and the 
number of operators needed depends strongly on the choice of the set of operators. How to choose these operators is 
a longstanding problem in nuclear physics E| . 

We select a small number of one-body operators F'^^\ k = 1, . . . ,n, assuming that the RPA eigenvectors can be 
approximated as linear combinations of the The approximate RPA vector is then given by 

n 

(27) 

k=l 

where J-'^^^ is the expectation value of F^'^'^ To determine the coefficients Ck the RPA matrices are projected onto 
the subspace {T^a^}: 

y^ki ^ jr(fc)5"/3^^^57'5jr(j)^ (28) 
^ki ^ jr(fc)B"/3jrW. (29) 
The RPA equation can then be expressed as 

M'^'q = (rin)^N'^'ci, (30) 

where h£L is a eigenfrequency of the projected RPA. The rank of the matrix we need to diagonalise to solve the 
RPA problem has been reduced from the number of 2-quasi-particle degrees of freedom to the number of one-body 
operators chosen. 
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D. Schrodinger equation on the collective path 



After having made a semi-classical approximation, which leads to a classical Hamiltonian, we need to remember that 
we are studying a quantum system. The standard technique to deal with this is to treat the classical Hamiltonian 
as a quantum one, and to calculate the eigenfunctions and energies. This is superficially similar to the generator 
coordinate method, especially in the Gaussian overlap approximation but it is actually rather different. The key 
point is the appearance of the kinetic terms, which correspond to time-odd generator coordinates (usually not include 
in the GCM). 

As discussed in Ref. 0, we can include all manner of quantum corrections to the potential energy, especially if we 
are interested in absolute values of the energy eigenvalues. On the other hand, shape mixing — a spread of the wave 
function along the collective path — is rather insensitive to these quantum corrections. Therefore, we shall consider 
the Hamiltonian along the collective path without further quantum corrections. 

One must include the zero modes when quantising the Hamiltonian, since they describe rotational and other 
excitations. Quantisation of the Hamiltonian in a metric coordinate space turns the kinetic energy into a Laplace- 
Bcltrami operator (see, e.g., Ref. |ll|) in the relevant space, 

H{X) = -^ + V(X). (31) 

The collective Schrodinger equation can then be written as 

i7(X)*(X) = £;*(X). (32) 

In this paper we discuss calculations with one true collective coordinate, and a number of additional momenta for the 
zero modes (denoted as pf^)'. two or three angular momenta, depending on whether the state is axial or not, and 
two operators connected to a change of phase of the proton and neutron pairing gap, associated with pairing rotation. 
These latter quantise as and j-^^- 

The potential V is invariant under all the zero-modes, and only depends on the collective parameter, 

V{X) ^ V{Q). (33) 

The Laplace-Beltrami operator, Ag, with variable (but diagonal) mass matrix , where the zero- mode masses are given 
by Bi ~ Bi{Q), can then be written as 

where g = Bq Bi, and Bq is identical to B in Eq. H15|l . and thus equals 1. Below we shall write B^^^^p for the 
neutron and proton pairing- rotational masses. These are calculated as 

B^. = nZV^^nII, (35) 
= Niy^Nll, (36) 

and the rotational moments of inertia are defined in the usual way. 

Since the potential and the masses are independent of the zero-mode coordinates the wave-function can be 
separated into various pieces, 

vE'(g,0Ar,0p,f]) =5-i/4{/(Q)-i=e™^"-i=^''*"^MK(f^)* (37) 

V ^TT V ^TT 

where m and k are the quantum numbers for neutron and proton pairing rotation, and /, M, K are the usual rotator 
quantum numbers. We shall be looking at ground states (bandheads) only, and therefore we shall now use I — M = 
K — 0, and since pairing rotational excitation corresponds to a change in particle number, we shall be use m = k = Q 
as well. Equation H34I) acting on 5* can now be rewritten as 



U"{Q) + ^""^ t''-''" c/(Q). (38) 
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Here we see the typical reason to absorb the factor g^l'^ into the wave function: it removes the linear derivative term, 
and we obtain a "centrifugal" potential in its place. This is fully consistent with the standard procedure for separation 
of variables in radially symmetric problems in two and three dimensions as can be found in any quantum mechanics 
textbook. Using (|37|l to separate variables, the Schrodinger equation can be written as 

Since we wish the wave function ^ to be normalisable we require it to be finite, and we must then insist that U{Q) 
goes to zero when g does. In the present work that only occurs when either of the pairing gaps collapses and thus 
B^p „ = is zero, and we shall ignore the rotational moments of inertia, which do not change very quickly. Below 
we solve Eq. H39(l on a grid with the boundary condition that [/(Qmax) = t^(Qmin) = 0. At points where i?^ = 
the condition U — holds exactly; for other cases applying this boundary condition will only give an upper limit on 
energy. 

The scaling of the wave- function removes g from expectation values, and the expectation value of any local operator 
A{Q) can evaluated as 

(.4) = J U{Q)A{Q)U{Q)dQ (40) 
which shows that U must be normalised according to 

j U{QfdQ = 1. (41) 



III. RESULTS 



To test the projection basis discussed in section FlI CI we implement our method for a interaction and configuration 
space that where the approximation can be compared with exact results. 



A. Pairing+quadrupole model 

We apply the LHA to the pairing+quadrupole Hamiltonian as described in . With a constraint on both neutron 
and proton numbers the Hamiltonian can be written as 

H' = H - ^IrNr, (42) 

T—n,p 

H = ^e.4cfc- J2 ^{PlPr + PrH)-^ jz QImQ^m, (43) 

k T=n,p A/=-2 

where tk are spherical single-particle energies, N-r is the particle number operator, Q2M is the dimensionless quadrupole 
operator 

Q2Af = -^E(^l^'^2A/|0 4c^ (44) 
V^Oq kl 

where &o = is the standard oscillator length and is the (dimensionless) pairing operator 

Pl-T.'y^k- (45) 

/c>0 

This Hamiltonian is treated in the Hartree-Bogoliubov approximation and it has been shown that at the minimum 
the local RPA for this Hamiltonian is equivalent to the quasi-particle RPA. For Af = 1, 2 we rewrite the quadrupole 
operators of Eq. H44|l as sums and differences, 

Q2M = ^ ± g2-M) , (M = l,2) , (46) 
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TABLE I: The quadrupole- and pairing interaction strengths, k and Gt, used for the examples discussed in this section. The 
deformation and pairing gap calculated for those interaction strengths arc also listed. 





K [MeV] 


G„ [MeV] 


Gp [MeV] 


e 


7 


A„ [MeV] 


Ap [MeV] 


^*Cr 


0.201367 


0.525586 


0.485390 


0.167 


0.0 


1.60 


1.60 


s^Fe 


0.122188 


0.379609 


0.478824 


0.213 


0.0 


1.60 


1.60 


e^Zn 


0.113687 


0.375498 


0.411877 


0.192 


0.0 


1.63 


1.70 


^''Zn 


0.113687 


0.375498 


0.411877 


0.286 


60.0 


1.95 


1.74 


^»Zn 


0.113687 


0.375498 


0.411877 


0.772 


0.0 


0.63 


0.96 



and the pairing operator of Eq. (|45|l as 

P^±) = (F, ± Pt) , r^p,n. (47) 

The pairing and quadrupole operators can then be arranged into five Hermitian, and four anti-Hermitian, Sj, 
operators: 

= {pi^\Pl,'-\Q2o,Q^fi\QiV), (48) 
s, = {Pt\Pt^,Qit\Qi2^)- (49) 
The Hamiltonian of Eq. (|43|l can then be written as 

H = Y^ €kclck -^Yl i^^R^Ri + \Y1 '^jSjSj, (50) 

k i j 

with — Gr for Pr^^ and — k for the Q operators. After solving the mean-field problem within the 

Hartree-Bogoliubov approximation the mass matrix and RPA potential around the minimum can be calculated as 

B"^ = Sa<5a0-2^Kj(S'j)a(5'j)/3, (51) 
j 

Va0 = S„<5„/3-2^K,(i?,)a(i?»)/3, (52) 
i 

where Ea = + Cqi is the 2 quasi-particle energy and a and /3 label 2 quasi-particle states. Oa is the 2 quasi-particle 
matrix element of the operator O, which can also be written as O^^/ with a = qq' . 

The spherical single particle energies are taken from 0- Our model space consists of two major shells. We follow 
and multiply all quadrupole matrix elements with the factor 



Nh 



C = Tr-^- (53) 

2 



where Nj^ is the harmonic oscillator quantum number of the lower major shell and Nh that of the higher one. To 
achieve the same root-mean-square radii for protons and neutrons different harmonic oscillator frequencies are adopted 
for each type of nucleons. As a result the proton and neutron quadrupole operators are multiplied by the factors 




and ap — \j (54) 

where N (Z) is the neutron (proton) number and A = N + Z. 

We have chosen a set of representative isotopes for the examples shown in this paper. For ^''Cr, ^^Fe and ^^Zn the 
interaction strengths are chosen to reproduce the ground state deformation listed in and the pairing strengths 
are chosen to approximately reproduce the relation give in The interaction strengths for the isotopes ^^Zn and 
^°Zn are chosen to be the same as for ^^Zn. The quadrupole- and pairing-strengths, k and Gt, are listed in table HI 
together with the corresponding deformations and pairing gaps. 
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TABLE 11: Comparing the full RPA energy, hQ, the projected RPA energy, hO,, and Sb.i for the (3-, 7-, A'^'- and A'^'-vibrations 
using the projection basis of 3- The energies are in units of MeV. 



/3- vibration 







hn 




5b 


Si 


54 Cr 


1.216 


1.246 




0.0013 


0.6528 


sspe 


2.511 


2.656 




0.0282 


0.0910 


s^Zn 


1.966 


2.033 




0.0053 


0.3593 


««Zn 


2.117 


2.188 




0.0050 


0.3117 


^»Zn 


1.029 


1.068 




0.0014 


0.2146 








7-vibration 








hfi 






5b 


5i 


54Cr 


2.289 


2.386 




0.0103 


0.0085 


58Fe 


1.971 


2.049 




0.0090 


0.0076 


^^Zn 


1.284 


1.298 




0.0010 


0.0010 


"^Zn 


1.280 


1.289 




0.0003 


0.0002 




3.166 


3.348 




0.0563 


0.0511 








A'^-'-vibration 








nn 


an 




5b 


Si 


54Cr 


3.208 


3.894 




0.1591 


0.1496 




3.239 


3.935 




0.2712 


0.2534 


«2Zn 


3.383 


4.173 




0.5220 


0.5107 


^'^Zn 


3.661 


5.161 




0.5171 


0.4618 


■^'^Zn 


1.814 


2.052 




0.0085 


0.0105 








A'^'-vibration 








nn 


hn 




5b 


Si 


54 Cr 


3.449 


4.505 




0.5237 


0.5113 




3.549 


4.824 




0.6442 


0.6060 


62Zn 


3.537 


4.797 




0.7934 


0.7751 


«6Zn 


4.255 


5.493 




0.7827 


0.7977 


^»Zn 


3.573 


4.824 




0.6801 


0.6809 



B. Improved approximate representation of the normal- mode operators 

The quality of the results achieved by the projection method described in section III CI strongly depends on the 
choice of the single particle operator basis. In Ref. it was demonstrated that a basis set consisting of pairing, 
multipole and spin dependent one-body operators are not able to reproduce the results of a full RPA calculation. On 
the other hand if the basis is chosen to be a set of state-dependent Hermitian one-body operators of the structure 

^^-E^(«^«^).+h-c. (55) 

a " 

good agreement can be achieved with a small set of operators. The suppression factor E^^ can be understood if 
one looks at a simple example . With the basis of Eq. H55|) good results can be achieved for the low- lying (3- and 
7- vibrations 0, as can be seen in table Ull with a small set of operators consisting of the 8 pairing and quadrupole 
operators 

Fk^ (^P^+\P^-\Q20r,Qitl) , T = n,p. (56) 

Even though the /3- and 7-vibrations are well described with this basis set, the higher lying solutions of pairing- 
vibrational character are not well described. A couple of examples are listed in table |n] These results are not 
significantly improved by including higher-order multipole or quadrupole-pairing operators in the basis. 

To improve the results for the pairing vibrations we include a pairing operator only active close to the Fermi-surface. 
To avoid the problem of having to select by hand which levels that would have a non-zero matrix element we simply 
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TABLE III: Comparing the full RPA energy, HQ,, the projected RPA energy, hD,, and Sb.i for the /?-, 7-, A*-^^- and A'^'-vibrations 
using the new projection basis 1571 . The energies are in units of MeV. 



/3- vibration 







hn 




5b 


Si 


54 Cr 


1.216 


1.225 




0.0003 


0.6633 


sspe 


2.511 


2.537 




0.0030 


0.0630 


s^Zn 


1.966 


2.004 




0.0023 


0.3613 


««Zn 


2.116 


2.131 




0.0010 


0.3156 


^»Zn 


1.029 


1.031 




0.0001 


0.2110 








7- vibration 








nn 






Sb 


Si 


54 Cr 


2.289 


2.386 




0.0103 


0.0085 




1.971 


2.049 




0.0090 


0.0076 


^^Zn 


1.284 


1.298 




0.0010 


0.0010 


^f^Zn 


1.280 


1.289 




0.0003 


0.0002 


^"Zn 


3.166 


3.348 




0.0363 


0.0311 








A'^-'-vibration 








nn 


hn 




Sb 


Si 


54Cr 


3.208 


3.212 




0.0008 


0.0022 




3.239 


3.247 




0.0021 


0.0039 


«2Zn 


3.383 


3.392 




0.0055 


0.0099 


^'^Zn 


3.661 


3.664 




0.0005 


0.0037 


■^'^Zn 


1.814 


1.815 




0.0001 


0.0088 








A'^'-vibration 








nn 


an 




Sb 


Si 


54 Cr 


3.449 


3.596 




0.0426 


0.0370 


5«Fe 


3.549 


3.647 




0.0176 


0.0309 


62Zn 


3.537 


3.621 




0.0210 


0.0201 


«6Zn 


4.254 


4.767 




0.9355 


0.9286 


^»Zn 


3.573 


4.039 




0.1649 


0.1855 



divide the standard pairing operator P± with a large power of Ea- If the suppression factor, E^, is chosen with a 
large enough k all matrix elements except the ones with Ea close to zero will become negligible and the result will 
not depend on k. The basis set is now 

Fk=iP^+\Pt\Q20r,Qitl^,^j, r^n,p. (57) 

We have chosen k — 10. From table IIIll we can see that almost all the low-lying vibrational modes are now described 
with a very high accuracy. 

To check if the wave functions are described as well as the RPA energies by the projection method we also calculate 
the overlap of the full RPA vector / and the projected RPA vector /. As criteria for good projection we use smallness 
of the following quantities 



Sb = 1- /,aS"^/"/3, (58) 

jfJl 

^JifJ) {ID 



S^ = l , (59) 



where 



a 



(60) 
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If ^ = the projection corresponds to an exact result. The difference between 5b and 5i is that an admixture of 
a spurious solution will contribute to 5i and not 5b] 5b is the consistent quantity from the topological analysis. In 
tableElwe can see that 5 has a relative small value for the /?- and 7-vibration but, as expected, a substantially larger 
value for the pairing- vibrations. The new projection basis does systematically improve the wave-function as well as 
the energy, as can be seen in table Hill where the values of 5 are much smaller. The exception being the second pairing 
vibration in ^^Zn. This is due to that the ordering of the RPA solutions in the full RPA relative to the projected RPA 
is different in this case. Since the improved basis for the projected RPA gives energy spectra and wave-functions that 
are much better than the set used in Ref. Q, we will use the new set in the following calculations. 

Calculating the collective path using the projection basis has an advantage besides reducing the rank of the RPA 
matrix. We do not have any spurious solutions in the projected RPA calculations since we have not included angular 
momentum and particle number operators in the basis. We can therefore avoid problems due to crossings between the 
spurious modes and the physical modes along the collective path. Away from the minimum there is still an admixture 
of the spurious solution into the collective coordinate. [As stated above, we should use a symplectic RPA to resolve 
this problem, but will not do so here due to its complexity.] 

C. Representative case for large amplitude collective motion 

We would like to perform as simple a test of our method as possible, especially, we would like the model space to be 
small. We decided to concentrate on ^^Fe; since this nucleus is 7-soft, it should provide a demanding testing ground 
for our methodology. We have calculated the collective path using both the full RPA and the projected RPA. 

The results are shown in Fig. HEl There is good agreement between the projected and full RPA results along the 
collective path in all cases. This provides a further confirmation of the quality of our projection basis and shows that 
the basis works well, also away from the mean-field minimum. The collective path is found in a smaller range of the 
collective coordinate when we are using the full RPA compared to the results for the projected RPA. This is due to 
the fact that the collective coordinate mixes with our spurious modes, which leads to problem with convergence in our 
double iterative method, due to the approximations made in the derivation. The mixing of the collective coordinate 
and the spurious solution remains small as long as the spurious mode is almost orthogonal to the collective solution. 
When the energy of the spurious solution is similar to the collective solution the denominator in the expression for the 
overlap becomes small which causes the overlap to become large. In the projected RPA calculation we do not have 
any spurious solutions since we have not included any of the operators connected with the spurious motion in our 
basis. Therefore we do not get a large spurious contribution to our collective coordinate and we have better numerical 
stability of our calculation. 

We first investigate axial collective motion (see also ^3)> by following the vibration (even though this is not the 
lowest eigenvalue at equilibrium, it is the lowest one of axial symmetry). From Fig. ^ we can see that the quadrupole 
moment is approximately proportional to the collective coordinate Q in the region —2 < Q < 0, which is an indication 
that we have a path relative close to what we would obtain from a mean-field calculation with a constraint on the 
quadrupole moment. At larger and smaller values of Q the deformation (Qo) remains almost constant. Instead, the 
collective coordinate is now dependent on the pairing fields, for large Q proton pairing and for small Q neutron pairing. 
At Q « 1.1 the proton pair field collapses to zero. Our collective path ends at this point, since the singularity at 
zero pairing is similar to the origin in polar coordinates, with Q playing the role of radial coordinate and the pairing 
phase (j) role of polar angle. The change from quadrupole to pairing mode is dominated by a narrowly avoided 
crossing with the lowest pairing-vibration at Q « 0.2. After this crossing the quadrupole moment, (Qo) saturates and 
the (Pp) starts changing. This avoided crossing shows that more then one collective coordinate would be needed for 
an accurate description of the collective dynamics. The RPA frequency of the /3-vibration is as expected proportional 
to the derivative of the cranking parameter, A. 

We have also looked at the potential energy, simply calculated as the expectation value of the Hamiltonian at each 
point. In Fig.^a) we see that the potential has a local energy maximum at Q ~ — 1, which corresponds to a spherical 
shape, and a shallow oblate minimum at Q ~ —1.6. The potential around the minimum show an a quadratic behaviour 
which indicates that the harmonic approximation in RPA is well fulfilled for small-amplitude collective motion, but 
obviously fails for wave functions that have substantial support away from the minimum. It can easily be seen that 
the regions where the potential energy has a positive derivative are the regions where the cranking parameter. A, has 
negative value and the converse. 

The key to the whole approach is the decoupling parameter, Z?, which is plotted in Fig. ^f). It has a small value 
indicating a good decoupling of the collective mode from all the non-collective modes. The exceptions are at Q w 0.2 
and Q < —2 which is due to two avoided crossings of the /^-vibration with the pairing vibrations, as can be seen in 
Fig-db). The large peak in Fig.^f) at Q > 0.5 is due to an approximate numerical over-completeness in in the basis 
on which we have projected the RPA vectors. The over-completeness comes when a pair field is zero and the projected 
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FIG. 1: (Color online) Large amplitude collective motion in **Fe with axial symmetry (following the /3-vibration, the second 
lowest RPA solution at equilibrium), a) Energy along the collective path, b) The square of the lowest RPA frequencies, c) 
The dimensionless quadrupole moment, (Qo), and the hexadecapole moment, (r'*y4o). d) The dimensionless pairing operators 
(Pt). e) The cranking parameter A and the chemical potential /i. f) The decoupling measure, D. The dotted line represents 
the numerically over-complete decoupling measure and the solid line the correct calculation where the contribution due to 
over-completeness of the projection basis has been removed. The grey (green online) curves represents the results for the 
projected RPA and the black curves are for the full RPA. 

mass matrix has a zero eigenvalue. This is due to the extra pairing term included in the basis in section Fill Bl The 
over-completeness appear near the collapse of the proton pairing field, as can be seen in Fig. ^d). Even tough the 
basis only becomes exactly over-complete at the point of pairing collapse the calculation of D is already influenced 
when the pair field is small, due to the fact the D is calculated from an inverse of the mass matrix (see Eq. (|24|) ') 
which becomes ill defined. We should of course remove such a spurious contribution; this can be done quite easily, 
and leads to the result plotted as a solid curve in Fig. ^f)- From now on we will only plot the value of D where 
the contribution from over-completeness has been removed. The collapse of the pair field has a surprisingly strong 
infiuence on the collective path. Whether this is a result of the approximations we made, our choice of force or a 
general feature is not clear at this point. 

In section Hi Dl it was described how to solve the one-body Schrodinger equation for the collective path. In the case 
discussed above the proton pair field collapses at Q ~ 1.1. We can therefore expect that proton pairing rotation will 
play a key role for the excitation spectrum of our system. The proton pairing mass is plotted in Fig.|21 We can see that 
close to the collapse of the proton pairing B^p oc {Q — Qmax)^ which is what we expect when the collective coordinate 
is approximately (Pp)- For negative Q B^p has a non-trivial behaviour. In Fig.|3we also show the lowest eigenvalue 
and radial wave function of the collective Hamiltonian including fc = proton pairing rotation. At Q = Qmax where 
B^p = the ground-state wave-function goes to zero linear in Q as expected from at the origin in polar coordinates. 
At Q = Qinin we have made the approximation that J7((3min) = 0. There is no bound state supported by the shallow 
oblate minimum and the lowest excited state is 1.30 MeV above the collective ground state. This excitation energy is 
substantially smaller then the RPA harmonic approximation energy of 2.54 MeV and reflects the anharmonic nature 
of the large-amplitude excitation. In table HVl we can see that the large component of the wave- function at small and 
negative Q gives rise to a reduction of the expectation value of Qq by almost 30% relative to the mean-field results. 
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FIG. 2: The upper panel show the proton rotational mass along the collective path for projected axial large amplitude collective 
motion in ^®Fe (following the second RPA solution). The lower panel show energy and the radial wave function for the large 
amplitude collective motion in ^**Fe, following the second RPA solution. The thick solid line is the potential energy, the thin 
horizontal line give the position of the lowest eigenvalue. The energy scale is on the left side. The wave function is shown as 
well (dotted line), with the scale on the right side. 



TABLE IV: The expectation values of the quadrupole- and pairing operators for k = collective Hamiltonian along the axial 
collective path in ^*Fe and at the mean-field minimum. 





Collective 


Klean — field 


(Oo) 


8.98 


12.37 


{Pn} 


4.36 


4.21 


(Pp) 


3.25 


3.34 



In the case where we follow the path emerging from the lowest mode, the 7-vibration, we can obtain similar results. 
A number of results identifying the collective path are shown in Fig.|21 where we have linear change of (Q2) with the 
collective coordinate, while all other expectation values remain relatively unchanged for \Q\ < 0.5. At larger values 
of the collective coordinate we see a saturation in (Q2) Siiid a strong reduction in the neutron pair field, which finally 
collapses to zero. Once again, this is mediated by an avoided crossing between quadrupole- and pairing-vibration 
modes. The decoupling measure, D, in Fig. |31f) has a similar behaviour as for the /3- vibration. The crossing with 
the pairing vibration is visible as an increase in D at around Q w 0.6. At larger Q we have a large contributions to 
D due to over-completeness of the basis this time caused by the strongly reduced neutron pairing field. By mirroring 
the potential to negative Q (and negative {Q2)) we get a closed collective path from the neutron pairing collapse at 
Q ~ 1.6 to the mirrored neutron pairing collapse at Q ~ — 1.6. The neutron rotational pairing mass in Fig. 0] shows 
the expected quadratic behaviour in \Q\ close to Qmax- In Fig. 01 we also show the eigenvalues and wave- functions of 
the collective Hamiltonian including to = neutron pairing rotation. The lowest excited state is at 1.49 MeV above 
the collective ground state. This excitation energy is substantially smaller then the corresponding RPA excitation 
energy of 2.05 MeV. This again is a result of the anharmonic nature of the collective potential. In tableElwe can see 
that the lowest state has a substantially reduced value of {Pn) compared to the mean-field value at the minimum. 



TABLE V: The expectation values of the quadrupole- and pairing operators for m = Q collective Hamiltonian along the 
non-axial collective path in ^^Fe and at the mean-field minimum. 





Collective 


Mean — field 


(Qo) 


12.63 


12.37 


(^^> 


3.53 


4.21 




3.31 


3.34 
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FIG. 3: (Color online) Large amplitude collective motion in ^^Fe (following the 7-vibration, the lowest RPA solution at 
equilibrium), a) Energy along the collective path, b) The square of the lowest RPA frequencies, c) The dimensionless 
quadrupole moments, (Qo), (Q2) and the hexadecapole moment, (r^y^o). d) The dimensionless pairing operators (Pt). e) The 
cranking parameter A and the chemical potential /j.. f) The decoupling measure, D. The grey (green online) curves are the 
results for the projected RPA and the black curves are for the full RPA. 




FIG. 4: The upper panel show the proton rotational mass along the collective path for projected large amplitude collective 
motion in ^^Fe (following the lowest RPA solution). The lower panel show energy and the radial wave function for the large 
amplitude collective motion in ^*Fe, following the lowest RPA solution. The thick solid line is the potential energy, the thin 
horizontal lines give the position of the lowest eigenvalue. The energy scale is on the left side. The wave function is shown as 
well (dotted line), with the scale on the right side. 
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TABLE VI; Comparison of the full RPA energy, fifi, the projected RPA energy, hQ, and Sb.i for the f3-, 7-, A'^'- and 
A*^ '-vibrations in case of i^«Gd and i«2os using the new projection basis 1571 . The energies are in units of MeV. 







HQ 




Sb 


Si 


/3-vibration 


0.8850 


0.9224 




0.0040 


0.5138 


7-vibration 


1.6860 


1.8490 




0.0366 


0.0325 


A'^'-vibration 


1.8089 


1.8147 




0.0049 


0.0042 


A^^'-vibration 


1.9923 


2.2256 




0.9409 


0.9339 








^«^0s 








nn 






Sb 


Si 


/3- vibration 


1.5704 


1.5851 




0.0058 


0.0087 


7-vibration 


1.1208 


1.1458 




0.0016 


0.0010 


A'^'-vibration 


1.6690 


1.6776 




0.0106 


0.0147 


A'^^'-vibration 


1.8476 


2.0673 




0.2165 


0.2826 



D. Realistic application of large amplitude collective motion 

The case of ^^Fe has the advantage that the configuration space is relatively small and therefor there are no big 
computational problems, and we could compare exact and approximate solutions. It still allows us to explore several 
key features of our method and test its feasibility and the quality of the results. To test the method in more realistic 
circumstances we decided to apply our method to the rare-earth region. We have chosen ^^^Gd and ^*^0s since the 
gadolinium nucleus is known to be /3-soft, whereas the osmium isotope is 7-soft. Both nuclei are situated in a region 
which is rich in nuclear structure phenomena. 

The calculations for ^'^^Gd and ^*^0s are done in a configuration space consisting of the N = 5(4) and 6(5) 
neutron(proton) shells. The spherical single particle energies and suppression factors of Eq. (|53|l and (|54|) arc again 
taken from Q. We compare the RPA energies and the RPA vectors calculated with the full RPA and using the 
projected approximation in table IVII We find a good agreement for the low-lying solutions. The second pairing 
vibration is somewhat too high in energy which is also reflected in a small overlap of the RPA vectors. The projection 
basis seems to work very well in the cases of heavier nuclei and larger configuration spaces examined here. 

In Fig. EHZI we can see the results of the large amplitude collective motion following the lowest axial symmetric 
solution. We have included both the results obtained with the full RPA as well those employing the RPA projected on 
a basis. In ^^^Gd the lowest solution is the /3 vibration but it also has quite large pairing components, as can be seen 
in Fig. [S] in the change in strength of the pair fields. There is in general a non-trivial structure of the collective path 
close to the mean-field minimum which also can be seen in Fig. |S1 a) where the energy show anharmonic behaviour 
around the minimum. Figure |Slf) show a good decoupling of the collective degrees of freedom from the non-collective 
degrees of freedom in the region close to the minimum. In Fig. we can see that there is one low-lying solution of 
the collective Hamiltonian. The fact the lowest eigenvalue is situated high in energy, relative to the range in which 
we have found the collective potential, tells us that the assumption that U (Q) = at the at the ends of the collective 
path is not justified in this case. One can expect that the correct U would stretch substantially outside the range on 
which we have calculated the collective path. 

^^^Os is a 7-soft nucleus and we show the result following the 2 lowest normal modes. In Fig.|7|we can see that the 
lowest axial RPA solution is mainly of proton pairing nature. The strength of the proton pair field is proportional to 
the collective coordinate and that the pair field collapses at Q « 1.6 which leads to a jump in the chemical potential. 
For small negative values of Q there is an avoided crossing with a mode that is dominantly a shape vibration, which 
leads to a reduction of (Qq). The energy along the collective path in Fig.|7|a) show a maximum when (Pp) and a 
approximately harmonic behaviour around the minimum. FigureElf) shows a good decoupling of the collective degrees 
of freedom from the non-collective degrees of freedom in the region close to the minimum. At large negative values 
of Q we have a crossing with a proton paring vibration which gives large state mixing and therefore no decoupling of 
the collective solution. 

Fig.|Slshow the results when following the collective path defined by the lowest 7-vibration in -'^^■^Os. The calculation 
shows that the collective path is mainly dominated by the increase of the (Q2) tri-axial deformation. At Q > 0.8 we 
see an avoided crossing of the /3- and 7-vibration which causes an numerical instability in our calculations. This also 
signals the need for more then one collective coordinate. Even though there are numerical difficulties in implementing 
our method in some cases we can see that our projection method works very well in the larger configuration spaces 
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FIG. 5: (Color online) Axial collective motion in ^^''Gd. See Fig. 0for more details. 




FIG. 6: The energy along the collective path for projected axial large amplitude collective motion in ^^®Gd (following the first 
RPA solution) is drawn with the thick solid line with the energy scale on the left side. The lowest eigenvalue of the Hamiltonian 
is drawn as a thin horizontal line. The corresponding eigenfunction has the scale on the right side. 



employed here and it is practically implementable. 



IV. CONCLUSIONS AND SUMMARY 

We have extended the method of calculating the self-consistent collective path presented in to include constraints 
on the particle number and implemented it for the quadrupole+pairing Hamiltonian 7] . The method consists of finding 
a series of points fulfilling the force equation, where the local direction of the collective path is determined in each 
point by the local normal modes. The local RPA equations and the force equation are solved in a double iterative 
process with constraints on the particle numbers and the step length along the collective path. The method allows us 
to determine the collective coordinate from the Hamiltonian without having to assume a priori which are the relevant 
degrees of freedom. 

To implement this method in heavier systems and for more realistic nuclear forces we need to truncate the RPA 
calculation in a way that will give an accurate approximation of the low- lying RPA solutions. We have improved the 
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projection method originally presented in p] in such a way that we are able to describe all low-lying states including 
pairing vibrations. This is done by expanding the state-dependent basis suggested in to include a pairing term, 
which is only active around the Fermi surface. The new basis set gives RPA energies very close to the exact solutions. 
The calculation of the overlap of the wave functions between the full and our approximate RPA solution shows that 
the wave functions are almost identical. We can therefore expect our method to give a good approximation to the 
collective path. 

Our method of calculating the collective path has been implemented for the cases of '^^Fe, ^^^Gd and ^^^Os. We 
have chosen to follow the lowest axially symmetric and tri-axial solutions. The decoupling of our collective coordinate 
from all the non-collective coordinates can be quantified in the decoupling measure D. This is found to be small 
along the collective path with the exception of regions of avoided crossings where the system undergoes configuration 
mixing. In such regions one collective coordinate is not sufficient to describe the system accurate. In regions where the 
projection basis gets over-complete D has to be calculated with special care. This happens when the proton and/or 
the neutron pairing collapses. 

We see that the collective path goes through avoided crossings with pairing solutions in most cases in both ^^Fe 
and ^^^Os. This leads to collapses of the pair fields and an end of the collective path. These avoided crossings also 
show that more then one collective coordinate would be needed for a accurate description of the collective motion. 
In the /9-soft nucleus ^^^Gd the collective coordinate is of a more complicated structure which can also be seen in 
the non-harmonic shape of the potential energy. Our projected local RPA method for calculating the collective path 
gives very good agreement with the results obtained using the full RPA. The method is very useful when calculating 
self-consistent large amplitude collective motion in large a configuration space. 

By solving the one-dimensional "radial" Schrodinger equation along the collective path we are able to examine the 
effect of the collective motion on ground state properties. In cases where the collective path ends with a collapse of the 
pair field we must include the effect of pairing rotations on the low energy spectrum. It is surprising that almost all our 
calculations are dominated by states with collapsing pairing; there may well be important lessons in this feature. The 
reason they occur so frequently is the presence of low- lying configurations without pairing. One might ask whether this 
is an artifact of our model, and whether larger configuration-spaces with more complicated interactions would behave 
differently. Such calculations are clearly called for, but we do not expect dramatically different results, since pairing 
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mainly acts in a small region around the Fermi-surface. There is also a slight possibility that the approximations we 
made in our treatment of spurious admixtures contributes to these effects. The surprising importance of the pairing 
collapse is the main result of our calculations that is usually not seen in a standard constrained mean-field calculation. 

In this paper we have implemented a method to find the adiabatic self-consistent collective path for a nuclei. A 
technique to truncate the basis in which the RPA equations are solved has been improved and a good agreement 
between the full and truncated RPA is found. To solve the RPA equations in a limited basis has proven to be a useful 
and practical way of calculating the collective path within the local harmonic approximation. We intend to apply 
similar techniques to the interesting problem of collective motion at finite rotational frequency in the near future. It 
remains to be investigated how we can include the covariant terms in the RPA in a suitable approximation. 
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Appendix 

The six conditions (Eqs. (|13|l and (|19H21|) ') gives a set of equations which can be solved for the six parameters e -"\ 
77jy', A^"' and /il""* for each i and n. For i — we get 

- ^ (Ql) 



1 



, , , , f(n) . ]v20(n) 
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For all other i Eqs. H16|l and (|17|l together with the constraints (|13H21() give 



(63) 



j-(n) )v20(") 
f(n) . f(n) _ { L 

_ gf(").Arg'(»)_A(")/(»).ivg'^(") 

^20(n) , ^2 



„(n) _ -^V^r - ^^ J' ' ' ff-f-^ 

f^^r - 20(n) ,.20(„) " 



These equations can easily be generalised to incorporate additional constraint operators like angular momentum. 
They are slightly more complicated than those shown in other work f^, since we have chosen to fix the step size in 
the collective coordinate. 
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